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Abstract 

A statistical analysis of the pre-SNO solar neutrino data is made assuming the solar 
neutrino deficit to be resolved by the interaction of the neutrino magnetic moment with 
the solar magnetic field. Given the general characteristics of the field profiles known to 
lead to the best event rate predictions, several specific choices of profiles are assumed 
and fits to the rates, to the electron recoil spectrum in SuperKamiokande and to the 
combined rates and spectrum are performed. Both the rate and spectrum fits are 
remarkably better than in a previous analysis of the magnetic moment solution which 
used a former standard solar model and previous data, and in particular, rate fits are 
far better than those from oscillations. The global analysis also reflects the general 
very good quality of the fits. 
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Whereas the apparent anticorrelation of the neutrino event rate with sunspot activity 
claimed long ago by the Homestake collaboration |IJ remained unconfirmed by other exper- 
iments @, @|, || and theoretical analyses IJ, the magnetic moment solution to the solar 
neutrino problem is at present an important possibility to be explored in the quest for an 
explanation of the solar neutrino deficit. This is the idea, originally proposed by Cisneros 
H, and later revived by Voloshin, Vysotsky and Okun 0, that a large magnetic moment 
of the neutrino may interact with the magnetic field of the sun, converting weakly active to 
sterile neutrinos. It now appears in fact that this deficit is energy dependent, in the sense 
that neutrinos of different energies are suppressed differently. In order to provide an energy 
dependent deficit, the conversion mechanism from active to nonactive neutrinos must be 
resonant, with the location of the critical density being determined by the neutrino energy. 
Thus was developed the idea of the resonant spin flavour precession (RSFP) proposed in 
1988 ||. It involves the simultaneous flip of both chirality and flavour consisting basically 
in the assumption that the neutrino conversion due to magnetic moment and magnetic field 
takes place through a resonance inside matter in much the same way as matter oscillations 
. A sunspot activity related event rate in a particular experiment would hence imply that 
most of the neutrinos with energies relevant to that experiment have their resonances in the 
sunspot range. However, the depth of sunspots is unknown (they may not extend deeper 
than a few hundred kilometers) and the observed field intensity is too small in sunspots to 
allow for a significant conversion. The anticorrelation argument has therefore lost its appeal 
for several years now. 

Despite the absence of the anticorrelation argument, several main reasons may be in- 
voked to motivate RSFP and investigating its consequences for solar neutrinos. In fact both 
RSFP and all oscillation scenarios indicate a drop in the survival probability from the low 
energy (pp) to the intermediate energy neutrino sector ( 7 Be, CNO, pep) and a subsequent 
moderate rise as the energy increases further into the 8 B sector. The magnetic field profiles 
providing good fits to the event rates from solar neutrino experiments typically show the 
characteristic of a sharp rise in intensity at some point in the solar interior, followed by a 
progressive moderate decrease , (Tl| . This is in opposite correspondence with the energy 



dependence of the probability in the sense that the strongest field intensities correspond to 
the smallest survival probabilities. Hence RSFP offers a unique explanation for the general 
shape of the probability, which naturally appears as a consequence of the field profile. On 
the other hand, from solar physics and helioseismology such a sharp rise and peak field 
intensity is expected to occur along the tachoclyne, the region extending from the upper 
radiative zone to the lower convective zone, where the gradient of the angular velocity of 
solar matter is different from zero [|T2 1 , [ 13| . Furthermore, it has become clear [14]], [ |10|| , 



TTj] that RSFP provides event rate fits from the solar neutrino experiments that are far 



better than all oscillation ones [15|, [16||. Finally, there are recent claims in the literature for 



evidence of a neutrino flux histogram containing two peaks, an indication of variability 
pointing towards a nonzero magnetic moment of the neutrino. 

The aim of this paper is to present a statistical analysis of seven field profiles proposed in 
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TCfl , all obeying the general features described above, in the light of the updated standard 
solar model (BP 2000) ||TH[ and the most recent data ]I|, 0, ||. For each profile the 
best fits of both the rate and recoil electron spectrum in neutrino electron scattering are 
determined along with the corresponding % 2 . Unlike a previous analysis, the contribution 
of hep neutrinos is taken into account in the spectrum fits. The results show a remarkable 



improvement in the quality of both sets of fits with respect to ref. |10|] where the BP '98 



standard solar model [19] and the available data at the time were used. Throughout the 
calculations a value of \i v = 10~ n jiR is assumed. 

The profiles investigated are the following: 

Profiles (1),(2) 

x <x R (1) 



B = 



B = B, 



x 



0" 



B = B n 



X — Xc 



Xr < X < X C 



Xc < x < 1 



1 - No- 
where x is the fraction of the solar radius and units are in Gauss. For profile 1 xr 
xc = 0.85 and for profile (2) x R = 0.65, xc = 0.80. 

Profile (3) 

B = , X < Xr 
X -Xr 



B = B, 



0" 



B = B 

with xr = 0.65, xc = 0.75. 
Profile (4) 



X C ~ X R 

x-0.7^ 2 



Xr < X < Xc 



0.3 



Xc < x < 1 



(2) 

(3) 
0.70, 

(4) 
(5) 
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B = B 
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x- 0.7^ 21 
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, X > Xr, 



with x R = 0.71. 
Profile (5) 



B 



B = , X < Xr 

Bo 



cosh 30 [X — Xr) 



,X>Xr 



with Xr = 0.71. 
Profile (6) 



B = 2.16 x 10 3 , x < 0.7105 



(7) 
(8) 

(9) 
(10) 

(11) 
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B = B 1 
B = 



:i p-0.75\ 2 ' 
~ V 0.04 J 

Bp 

cosh30(x - 0.7483) 



0.7105 < x < 0.7483 
, 0.7483 < x < 1 



with B Q = 0.998^. 
Profile (7) 
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B = Bn 



B = 2.16 x 10 3 , x < 0.7105 
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1 - 



0.04 



0.7105 < x < 0.7483 



B = 1.1494B [1 - 3.4412(a; - 0.71)] , 0.7483 < x < 1. 



(14) 
(15) 
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The ratios of the RSFP to the SSM event rates R% a 

ci sk are defined as before [|1(| with 
the exception of the SuperKamiokande one in which the energy resolution function |2(J is 
now taken into account 
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th 
SK 



Ei Jo" dE e j£« d£5/(^ e> £?.) /|f d£<M£)[P(£)^ + (1 " P(E))^ 

em 

E< /o°° dE e dE>f(E> ei E e ) f§£ dEHE)^- 



(17) 



Here (j>i{E) is the neutrino flux for component i (i = hep, 8 B) and f(E' e ,E e ) is the energy 
resolution function of the detector in terms of the physical (E e ) and the measured (E e ) 
electron energy (E e = T + m e ). The lower limit of E e is the detector threshold energy 
(E em = E eth with E eth = 5.5MeV) and the upper limit is evaluated from the maximum 
neutrino energy E M 



T 



2E 2 



M 



M 



m e + 2E 



M 



For the lower jKJ and upper |2T|] integration limits of the neutrino energy one has respectively 



E n 



T' + \T' 2 + 2m e T 



, E M = 15MeV (i = 8 B) , E M = 18.8MeV (i = hep). (19) 



The weak differential cross sections appearing in equation (17) are given by 



daw GF 2 m, 



dT 



_ [(g v + qa) 2 + (gv - 9a) 2 (1 - F 
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(20) 



for v e e scattering, with gy = \ + 2sin 2 9w, 9a = h- For v^e and v T e scattering, 



Iv -9a 
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E 2 J 



(21) 
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with g v = — \ + 2sin 2 9w, 9a = — §■ 

The recoil electron spectrum in SuperKamiokande is now defined as 



Rf K = a* -3 (22) 



with i = hep, 8 B for 18 energy bins (j=l,...,18) @. 

The fluxes and partial event rates for each neutrino component in each experiment 
were taken from |18| and the solar neutrino spectra from Bahcall's homepage f2~I|| . The 



contribution of the hep flux to both the Gallium and Homestake event rates was neglected. 
The x 2 analysis for the ratios of event rates and electron spectrum in SuperKamiokande 
was done following the standard procedure described in [|10J [[ 

The ratios of event rates to the SSM event rates and the recoil electron spectrum normal- 
ized to the SSM one, both denoted by R th in the following, were calculated in the parameter 
ranges Am^ = (4 — 22) x lO^eV 2 , B = (3 — 30) x 10 4 G for all magnetic field profiles and 
inserted in the \ 2 definitions for the rates and spectrum, 

Urates = E (i# ~ ~ RjS*) (23) 

iij2=i 1 2 

with (Ga=l, Cl=2, SK=3) 

18 

X% = E (R% ~ K P ) KM)A (*£ - RT) ■ (24) 
h j'2=i 

The quantities R exp in eqs. (21), (22) are directly read from tables I, II respectively and 
the total error matrices a 2 (tot) are derived from the definitions given in [^Oj] , using |21j and 



the error bars in tables I, II. In fitting the rates one has 3 experiments and 2 parameters 
(the mass square difference between neutrino flavours and the value of the field at the peak), 
hence the number of degrees of freedom is one, while for the spectrum there are 18 data 
points, hence 16 degrees of freedom. For global fits one has 

X% = E ( R n- WBNSiW - K P ) ■ (25) 
h J2=i 

Here index values 1, 2, 3 are used for the rates and 4, 21 for the spectrum. The best fits 
for the rates in terms of Am^Bo and the corresponding values of \ 2 f° r rates, spectrum 
and global fits are shown in table III. In table IV the spectrum best fits are shown with 



^Here only the main definitions and differences are registered. For the calculational details we refer the 



reader to ref. 10 
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the corresponding values of y 2 „„ and table V shows the global best fits with the values of 
Xratesy X% i n eacn profile. In this case the number of d.o.f. is 19 (=21-2). 

It is seen that the rate fits (table III) are in some cases excellent, with four of the profiles 
being much favoured, in particular profile 2. The most unfavoured ones (profiles 4, 5 and 7) 
give nevertheless substantially smaller x 2 as compared to oscillations |15], |IEfl . It is difficult 
at present to develop a general distinction between the characteristics of those profiles which 
provide the best fits and the others. In profiles 6 and 7 for instance the field intensity along 
the bottom of the convective zone rises in the same way and the quality of the rate fits is 
quite different: this suggests that an upward facing concavity (profile 6) along the upper 
convective zone is preferred relative to a constant slope (profile 7). In fact, since the low 
energy sector of the 8 B neutrinos is highly suppressed because their resonances are located 
close to the intermediate energy ones (pep, CNO) whose suppression has to be maximal, 
the higher energy sector of 8 B has to undergo little suppression in order to ensure a survival 
probability close to 1/2 as reported by SuperKamiokande. Hence the magnetic field should 
decrease steeply at first after its maximum and then more smoothly. The same conclusion 
is also suggested from inspection of profiles 4 and 5: while they exhibit a sudden rise at 
the bottom of the convective zone with infinite slope, their downward slopes have opposite 
shapes to each other, favouring the case with an upward facing concavity (profile 5). As 
for profile 4 the fits still worsen, as expected, if the exponent 2 in eq. (8) is increased. An 
infinite rise in field intensity (profiles 4 and 5) is also disfavoured relative to a finite one 
(compare 3 with 4 and 6 with 5). The case of the two triangle profiles (1 and 2) suggests 
that the upward slope should be higher than the downward one. 

Table III shows unusually large values of \ 2 g i at the rate best fits for profiles 4 and 7. This 
is related to a high instability of x 2 both for rates and global fits against small variations of 
Am^ and B not shared by other profiles, as can be seen from a comparison with table V. 

The 90 and 99% CL regions for the rates defined in the Am^ , B plane from the con- 
dition x 2 < Xmin + ^X 2 with, for one degree of freedom, A% 2 = 2.71 and 6.64 respectively, 
are shown in fig.l. The four profiles chosen are those for which Xmin < 0.1. As far as the 
spectrum fits are concerned, the results for the minima of x% are a l so remarkably good 
(see table IV). A comparison of tables III and IV and inspection of fig.l show a discrep- 
ancy between the rate and spectrum best fits. To further examine this discrepancy one can 
analise for each profile the difference between the CL's corresponding to x% &t the rate and 



spectrum best fits also listed in table IV p2| . For 16 d.o.f. it is seen from the values of this 
difference that all rate best fits lie within a region of 10.5% CL with respect to the spectrum 
best fit: the function x% has a weak dependence on the parameters B and Am^, so that 
the mismatch between the rates and spectrum fits is statistically meaningless. Hence one 
should expect very good global fits. In fact the global best fits, all with a remarkably low 
X 2 g iobai ( see table V) always lie within a similar (11.4% CL) region relative to the spectrum 
ones and also very close to the rate best fits, namely within a 41% CL region with respect 
to these, as can be derived from tables III, IV and V and using ref. p2f. 
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As an example, the predicted spectrum (eq. (22)) for profile 1 is shown in figs. 2 and 
3 respectively at the rate and spectrum best fits superimposed on the SuperKamiokande 
data. The moderate rise occuring for E e > YIMeV is the effect of hep neutrinos. 

As previously referred, a value of //„ = 10~ n /i,B wa s assumed in the calculations. Since 
the order parameter is the product fiB and the peak value of the field is probably as large as 
3 x 10 5 G |L2|, |T3|1 , the best fits imply an effective value of \i v in the range (2 — 5) x 10~ 12 fj, B 



in consistency with most astrophysical bounds 

To summarize, the analysis of prospects for the magnetic moment solution to the solar 
neutrino problem made on the basis of BP'2000 solar standard model and using the most 
recent pre-SNO data reveals extremely good fits for the rates for a class of field profiles with 
a steep rise across the bottom of the convective zone and a more moderate decrease up to 
the surface. There is a substantial and noteworthy improvement in the quality of the fits 
with respect to a previous analysis using BP'98 |IJ and previous data. Spectrum fits are 
also quite good and consistent with the rate fits. 

In the present analysis only time averaged data were considered and a fitting was made 
to a time constant profile 'buried' in the solar interior. If, on the contrary, the active 
neutrino flux turns out to be time dependent, a situation most likely to be interpreted 
through the magnetic moment solution with a time dependent interior field, the present 
approach is obviously inadequate and all fittings thus made will in some way be flawed. 
Averaging the event rates over time implies disregarding possible information in the data 



which otherwise is available if different periods of time are considered [17]. The robustness 



of such a procedure will greatly improve with the accumulation of more data. 
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Experiment 


Data 


Theory 


Data/Theory Reference 


Homestake 


2.56 ±0.16 ±0.15 


7.7±i;i 


0.332 ±0.05 § 


Ga 


74.7 ±5.13 


129±g 


0.59 ±0.06 ||,|§ 


SuperKamiokande 


2.4 ±0.085 


5.15±J;? 


0.465 ±0.052 § 



Table I - Data from the solar neutrino experiments. Units are SNU for Homestake and 
Gallium and I0 6 cm~ 2 s~ 1 for SuperKamiokande. The result for Gallium is the combined one 
from SAGE and Gallex±GNO. 



Energy bin (MeV) R £ a 



5.5 < E e < 6 


0.461 ±0.025 


6<E e < 6.5 


0.435 


!±0.02 


6.5 < E e < 7 


0.455 


±0 


.016 


7 <E e < 7.5 


0.474 


±0 


.016 


7.5 < E e < 8 


0.490 


±0 


.018 


8 < E e < 8.5 


0.477 


±0 


.018 


8.5 < E e < 9 


0.455 


±0 


.018 


9 < E e < 9.5 


0.463 


±0 


.021 


9.5 < E e < 10 


0.472 


±0 


.026 


10 < E e < 10.5 


0.460 


±0 


.027 


10.5 < E e < 11 


0.452 


±0 


.027 


11 < E e < 11.5 


0.473 


±0 


.033 


11.5 < E e < 12 


0.457 


±0 


.037 


12 < E e < 12.5 


0.429 


±0 


.041 


12.5 < E e < 13 


0.488 


±0 


.049 


13 < £ e < 13.5 


0.493 


±0 


.058 


13.5 < E e < 14 


0.583 


±0 


.062 


14 < £ e < 20 


0.505 


±0 


.078 



Table II - Recoil electron energy bins in SuperKamiokande (1117 days) and the corre- 
sponding ratio of the experimental to the SSM event rate 0. 
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Profile Best fit (Amj!, Bp) xL^/ld-o.f. xyi6d.o.f. x^/lgd-ofc 

1 (7.08 x 10- 9 e1/ 2 , 6.6 x 10 4 G) 

2 (1.25 x 10- 8 eV 2 , 1.28 x 10 5 G) 

3 (1.26 x l(T 8 ey 2 , 9.5 x 10 4 G) 

4 (1.7 x 10- 8 eV 2 , 1.7 x 10 5 G) 

5 (2.1 x lirW 2 , 1.42 x 10 5 G) 

6 (1.65 x !0- 8 eV 2 , 9.7 x 10 4 G) 

7 (1.50 x 10- 8 eF 2 , 2.3 x 10 5 G) 

Table III - Rate fits: the values of Ami and B at the best fit of the rate for each 



3.71 x 10" 


-2 


9.85 


12.73 


9.13 x 10" 


-4 


11.1 


12.82 


5.96 x 10- 


-2 


9.72 


13.24 


0.357 




11.45 


41.46 


0.133 




11.7 


11.93 


9.51 x 10- 


-3 


11.35 


14.19 


0.38 




10.94 


81.44 



b 21 

profile, the value of Xrates min an d the corresponding values of x% and Xtot- 



Profile Best fit (Am 2 , , B ) \; ; , /Uid.o.r. CL(x 2 sPm J- CL( X % rn1 



1 (3.7 x 10- 9 eV 2 , 4.3 x 10 4 G) " ~ 9.46 1.89% 

2 (6.4 x 10- 9 eF 2 , 1.26 x 10 5 G) 9.93 6.70% 

3 (7.7 x !0~ 9 eV 2 , 9.5 x 10 4 G) 9.48 1.15% 

4 (< 10- n eV 2 , 8.66 x 10 4 G) 9.74 9.90% 

5 (9 x 10- 9 eF 2 , 1.43 x 10 5 G) 9.94 10.5% 

6 (7.8 x lO- 9 eV 2 , 9.5 x 10 4 G) 10.0 7.93% 

7 (< lO- 1 ^^ 2 , 2.27 x 10 5 G) 9.75 6.62% 



Table IV - Spectrum fits: the values of Am^ and B at the best fit of the spectrum 
for each profile, the corresponding value of X% min an d the difference between the confidence 
levels (or goodness of fits) of the spectrum and the rate best fits for 16 d.o.f.(see the main 
text for details). For profiles 4 and 7 x% i s an asymptotic value for small Am^. 



Profile 


Best fit (A? 


■n 2 21 ,B ) 


xL,/19d-o.f. 


X 2 rates/ ld -0-i- 


X?p/16d.o.f. 


1 


(7.1 x 10- 9 eV 2 , 


6.7 x 10 4 G) 
1.27 x 10 5 G) 
9.4 x 10 4 G) 
1.69 x 10 5 G) 


9.94 


0.100 


9.79 


2 


(1.32 x 10- 8 e^ 2 , 


11.83 


0.290 


11.33 


3 


(1.25 x 10~ 8 eV^ 


9.81 


0.116 


9.63 


4 


(1.67 x !0- 8 eV 2 , 


12.7 


0.708 


10.95 


5 


(2.13 x 10- 8 eV 2 , 


1.44 x 10 5 G) 
9.7 x 10 4 G) 
2.28 x 10 5 G) 


11.88 


0.133 


11.82 


6 


(1.78 x 10- 8 eF 2 : 


12.29 


0.366 


11.83 


7 


(1.45 x 10~ 8 eV 2 , 


12.0 


0.764 


10.16 



Table V - Global fits: the values of Am^ and B at the global best fit for each profile, 
the corresponding values of xfates and xlpectmm- 
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Figure captions 

Fig. 1. 90% and 99% CL regions with respect to rate best fits on the Am^, B (peak 
field value) plane. Units are eV 2 (x-axis) and Gauss (y-axis). From left to right and top 
to bottom: profiles 1, 2, 3, 6. Circles denote the rate best fits and diamonds the spectrum 
best fits (for the corresponding values of x 2 see tables III and IV). The apparent mismatch 
between spectrum and rate best fits is statistically meaningless as can be seen from the last 
column of Table IV and the global fits (Table V). 

Fig. 2. The electron recoil spectrum prediction for profile 1 for the values of Am^, 
Bo corresponding to the rate best fit (table III) as a function of the total electron energy 
superimposed on the data set for 1117 days [0. For this curve x 2 P — 9.85 (see table III). 
Units are in eV. Notice the moderate rise in the theoretical curve for E e > YIMeV from 
the effect of hep neutrinos. 

Fig. 3. Same as fig.2 for the values of Am^, B corresponding to the spectrum best fit. 
For this curve x% = 9-46 (see table IV). 
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